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A UNIQUENESS CRITERION FOR VISCOUS LIMITS OF 
BOUNDARY RIEMANN PROBLEMS 

CLEOPATRA CHRISTOFOROU AND LAURA V. SPINOLO 

o 

Abstract. We deal with initial-boundary value problems for systems of con- 
servation laws in one space dimension and we focus on the boundary Riemann 
Cn , problem. It is known that, in general, different viscous approximations provide 

Cn ' different limits. In this paper, we establish sufficient conditions to conclude 

that two different approximations lead to the same limit. As an application 
of this result, we show that, under reasonable assumptions, the self-similar 
second-order approximation 

,S^ • and the classical viscous approximation 

provide the same limit as e — > 0"*". Our analysis applies to both the character- 
istic and the non characteristic case. We require neither genuine nonlinearity 
nor linear degeneracy of the characteristic fields. 
Date: July 23, 2010 

^T) ', 1. Introduction and main results 

j«^ . We are interested in the system of conservation laws 

l> ■ dtU + d,F{u) ^ 0, (1.1) 

f— ^ ' where the unknown function U{t,x) takes values in M". the variables t and x are 

both scalar and the flux F : R" — )• M" is a smooth function. We assume that 
system (jl.ip is strictly hyperbolic, namely for every U g R" the Jacobian matrix 
DF{U) has n real and distinct eigenvalues 

Ai([/)<A2([/)<---<A„([/). 

For an exposition of the current state of the theory of systems of conservation laws, 
we refer to the books by Dafermos [S] and Serre |28j . 

In this paper, we study initial boundary value problems and we establish condi- 
tions ensuring that two different viscous approximations lead to the same solution 
in the limit. Then, we apply our result to self-similar second-order approximations 
and the classical viscous approximations. 

To contextualize the problem, let us first consider the Cauchy problem obtained 
by coupling (|l.ip with the initial datum 

U{{),x)^Uc{x). (1.2) 

The so-called Riemann problem is posed in the case that Uc{x) has the form 

Uc(x) ~ < ,, .^ „ 1.3 

^ ' \ U li X < 0, ^ ' 
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where U^ , U^ are two given constant states in R". The Riemann problem ()1.1[) . ()1.3p 
has been extensively studied and its analysis provides important information on 
both the local (in spacc-tinic) and the long-time behavior of the solution of a gen- 
eral Cauchy problem (|l.ll) - (ll.2p . Also, it serves as a building block for the con- 
struction of different approximation schemes (e.g. the random choice method of 
Glimm [13], the wave front-tracking algorithm, see Bressan et al [6] and Holden 
and Risebro |15j ) that lead to existence and uniqueness results for general Cauchy 
problems. It should be noted that these results hold for general flux F under the 
hypothesis that the total variation of the Cauchy datum Uc{x) is small enough. 
Suitable counterexamples ensure that handling data with large or unbounded total 
variation requires the assumption of more restrictive conditions on the structure of 
the flux F in (jl.ip . In the present paper, we deal with general flux functions and, 
therefore, we focus on Riemann data C/"*", U~ in ()1.3|) that are sufficiently close. 

One of the major difficulties posed by the Riemann problem is the non- uniqueness 
of distributional solutions. Various selection principles, often motivated by physical 
considerations, have been introduced in an attempt at singling out a unique solu- 
tion: see, for example, the entropy admissibility criterion and the conditions named 
after Lax [22] and Liu [23l [25] . A solution of the Riemann problem (|1.1|) , p.3p sat- 
isfying suitable admissibility criteria was first constructed by Lax [52] , Liu [53] and 
later by Tzavaras |31j and Bianchini [3]. 

In particular, in [3], solutions of (|l.ip . (|1.3p are obtained by considering the 
viscous approximation 

dtU' + d.,F{U') = e d.,{B{W)d,W) , (1.4) 

where _B is an n x n matrix which depends on the physical model under consid- 
eration and e is a positive parameter. As e -^ 0+, the family of functions [/^ is 
expected to converge to a distributional solution of (jl.ip . It is a challenging open 
problem to establish a rigorous proof of the convergence U^ —?' U in the general 
case, but results have been achieved in specific cases. See Bianchini and Bressan [i] 
and the references in the books [9] [28] . It should also be noted that the analysis 
in [3] ensures that the function constructed in [3] |22l [24] [31] is the unique solu- 
tion of (jl.ip . ()1.3p which is self-similar and can be obtained by patching together 
a countable number of rarefaction waves and of shocks (or contact discontinuities) 
satisfying the admissibility condition named after Liu. As a consequence, the limit 
as e — > 0+ of the family of functions C/^ in (|1.4p does not depend on the choice of 
B. As we see in the following, this is not true in the case of initial-boundary value 
problems. 

Another approach to the analysis of the Riemann problem was introduced inde- 
pendently by Dafcrmos [S], Kalasnikov [T^ and Tupciev [5S] and is based on the 
analysis of the family of parabolic problems 

dtU' + d^F{W) = etd:,{B{W)d:,W), (1.5) 

coupled with the initial datum 

Because of the presence of the "t" factor in the second order term, the Cauchy prob- 
lem (|1.5|) - ()1.6|) admits self-similar solutions in the form U'^{t,x) = V'^{x/t). Com- 
pactness results for the solutions of ()1.5|) - ()1.6|) have been established under suitable 



UNIQUENESS CRITERION FOR THE BOUNDARY RIEMANN PROBLEM 



hypotheses on B. See in particular Tzavaras [30l[3T], Joseph and LeFloch [18] and 
the references in Dafermos [H Section 9.8]. 

In the present paper, we consider the initial-boundary value problem for the 
system of conservation laws (|1.1|) . hence we assume that cc > 0. The analogous 
to the Riemann problem is the so-called boundary Riemann problem, obtained by 
coupling (jl.ip with the Cauchy and Dirichlet data 

U{0,x):^Uo, forx>0, C/(t,0)==C7, for i > 0, (1.7) 

where Uq and U are two given constant states. Since we are interested in solutions 
of small total variation, we restrict to the case that \U — Uo\ is small. The analysis 
of the boundary Riemann problem provides again information on the local and the 
long-time behavior of the solution of a general initial-boundary value problem and, 
moreover, it constitutes the building block for the construction of approximation 
schemes that lead to general existence and uniqueness results. In addition, because 
of the presence of the boundary, additional challenges are raised compared to the 
Riemann problem. First of all, the initial-boundary value problem (jl.ip . (|1.7p is in 
general ill-posed, i.e. it posseses no solution. Therefore, a notion of admissible set 
of boundary values can be introduced, see Dubois and LeFloch [TU] . 
Other challenges arise when studying the viscous approximation 

dtU' + d^F{U^) = e d^{B(U')d^U') 

C/^(0,a-) = L/o, x>0 (1.8) 

C/^(t,0) = L/d, t>0. 

It should be noted that, if the matrix B is invcrtible, the above initial-boundary 
value problem is well-posed for every Udi while if B is singular, one has to intro- 
duce a more complicated formulation of the boundary condition, see for example 
Bianchini and Spinolo [5] . We emphasize that B is singular in most of the physi- 
cally relevant examples, but to simplify the exposition, we focus on the case of B 
invcrtible. 

The first difficulty one has to address when studying (|1.8|) is the following: as- 
sume that as e ^ 0+ the family [/^ converges, in a suitable topology, to some limit 
U . Also, assume that U is self-similar function, [/(i, x) = V{x/t) and that V has 
bounded total variation, so that lim^_j.o+ V^(C) is well defined. Because of boundary 
layer phenomena, in general 

f/(t,0)= lim V{O^Ud. 

Also, as pointed out by Gisclon and Serre [H] [12], in general, if F, Uq and Ud 
in (|1.8p arc fixed and the matrix B varies, then the limit U varies. To understand 
this behavior, let us first focus on the case that the boundary is non characteristic, 
namely all the eigenvalues of the Jacobian matrix DF(U) are bounded away from 
and let us denote the trace by 

[7= lim V(£,). 

See condition (|1.20p in Section 11.11 for the rigorous definition of non characteristic 
boundary. Then one expects that there is a boundary layer of (|1.4[) connecting U 
with Ud, namely the system 

B{W)W' = FiW) - F{U) 

W{0) = Ud \m-iy^+oc W{y) = U ^''^^ 
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admits a solution W{y). Here, ' denotes the first derivative. Since the solvabil- 
ity of (|1.9p depends on B, then also U does. The boundary characteristic case 
occurs if an eigenvalue of the Jacobian matrix can attain the value 0, see condi- 
tions (jl.2ip - p.22p in Section [01 for the precise definition. The boundary case is 
more complicated to handle because of the possible interactions between boundary 
layers, traveling waves and contact discontinuities, see the analysis in Ancona and 
Bianchini [2] and Bianchini and Spinolo 0. However, one again obtains that in 
general the limit depends on the viscosity matrix B. 

In |17j . Joseph and LeFloch studied the viscous approximation 

dtU' + d,FiU') = etd,{B{U')d,U') 

U''{0,x) = Uo. x>0, (1.10) 

U%t,0)^UD, t>o, 

in the case that B is the identity matrix, B{U) = I. They established compactness 
results and provided a detailed description of the limit in both the case of a non 
characteristic and (under some technical assumptions) a characteristic boundary. 
In [TS], among other results, Joseph and LeFloch extended the analysis to more 
general viscosity matrices. 

The motivation of this work is to investigate whether the limit of the classical 
viscous approximation (jl.Sp and of the self-similar viscous approximation (|1.10|) 
coincide. In the case of the Riemann problem, the limit of ()1.4p . (|1.6p and of 1)1.5^ - 
(jl.6|) are expected to coincide, because of the uniqueness result in Bianchini [3]. 
However, the case of the boundary Riemann problem is a priori much less clear, 
because of the work by Gisclon and Serre [TTJ [H] stating that, in general, the 
limit of the viscous approximation depends on the viscosity matrix. On the other 
hand, the analysis in Joseph and LeFloch [T71 [TH] ensures that the equation of the 
boundary layers of the self-similar approximation (|1.5p is the same as the equation 
of the boundary layers of the classical viscous approximation (|1.4p . namely it is 

[BiU)U']'^[FiU)]' . (1.11) 

Hence, to determine whether the limits of (jl.Sp and (|1.10p coincide, a first step is 
to investigate whether there exists a unique solution of the limiting system which 
can be constructed by patching together boundary layers satisfying (|l.lip and a 
countable number of rarefaction waves and shocks (or contact discontinuities) that 
arc admissible in the sense of Liu. 

Our first result provides an affirmative answer to the previous question. Here, 
we state the main result of this paper: 

Theorem 1.1. Assume that Hypotheses 1-3 given in Section \l.l\ hold, let Uo,Ud 
be two given constant states in R" and consider the system of conservation laws 

dtU + d^F{U)^0. (1.12) 

Then, there exist positive constants 6 and C , S sufficiently small, such that, if \Uo — 
Ud\ < (5; then there exists a unique distributional solution U of (|1.12p satisfying 
the following properties: 

(1) the function U attains the Cauchy datum U{0,x) ~ Uq for a.e. x; 

(2) the function U is self-similar, namely U{t,x) — V{x/t) for a suitable func- 
tion V which also satisfies TotVar V < CS; 

(3) all the shocks and the contact discontinuities of U are admissible in the 
sense of Liu; 



UNIQUENESS CRITERION FOR THE BOUNDARY RIEMANN PROBLEM 5 

(4) let I be any open interval where V is continuous, then V is differentiahle 
in I; 

(5) there exists a value U_ with \Uq ^ U\ < CS, such that the following two 
conditions hold. 

(a) // we denote by U the trace of U on the t-axis, namely 

U := lim V(^), (1.13) 

then F{U_) ~ F{U) and the shock or the contact discontinuity connect- 
ing U_ (on the left) with U (on the right) is admissible in the sense of 
Liu. 

(b) There exists a boundary layer W satisfying 

B{W)W' = F{W) ^ F{U) 

W{0) = Ud limy^+oo W{y) ^U ^ ' ' 

and \W{y) ~ U\ < CS, \W'iy)\ < CS for every y. 

It should be mentioned that Liu's admissibihty condition is stated in Section r2.31 
see formula (j2.14p . For an extended discussion, see also Dafermos [SI Chapter 8.4]. 

Some remarks are here in order. First, the novelty in Theorem 1 1.1 1 is the unique- 
ness part, since the existence part is already known. Second, if the boundary is 
non-characteristic then F{]f) — F{U) implies that U_ = U, hence (|1.14p reduces 
to (|1.9p and the proof of Theorem 1 1 . 1 1 is easier. In the boundary characteristic case, 
one has to handle the possibility of a zero-speed shock (or contact discontinuity), 
so the value U_ comes into play and the analysis is considerably more complicated. 
Last, the hypotheses introduced in Section [TTT] are quite standard in this context 
and they imply that the viscosity matrix B is strictly stable in the sense of Majda 
and Pego |5B] (see Corollary 2.2 therein) and they are the same assumptions con- 
sidered by Gisclon and Serre in [11] [12]. In particular, we require neither genuine 
nonlincarity nor linear degeneracy of the characteristic fields. Our hypotheses im- 
ply that the viscosity matrix B is invertible. The extension to the case of a singular 
viscosity matrix, which is more interesting in view of physical applications, raises 
no apparent difficulty provided that the assumptions introduced by Kawashima 
and Shizuta }21) and a condition of so-called block linear degeneracy introduced 
in Bianchini and Spinolo [5] are all satisfied. However, in the present paper, to 
simplify the exposition, we restrict to the case of an invertible viscosity matrix. 

An application of Theorem 11.11 to the classical and the self-similar viscous ap- 
proximation shows that the two limits coincide even in the case of initial-boundary 
value problems. Here is the result: 



Theorem 1.2. Let Hypotheses 1-3 in Section \l.l\ hold, then there exist constants 
C and 5, S sufficiently small, such that, given two values Uq, Ud G R" satisfying 
\Uo — Ud\ < S, the following holds. Consider the classical viscous approximation 

dtU' + d,F{U') = e d^{B{U')d.,U') 

U%0,x) = Uo, x>0, (1.15) 

U^t,0)^UD, t>0, 
and the "self-similar" approximation 

dtZ^ + d^F{Z^) = e t d,{B{Z^)d^Z^) 

Z^(0,a;) = [/o, a;>0, (1.16) 

Z^{t,Q) = UD, t>0, 
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and assume that as e ^ 0+ 

W^U and Z^ ^ Z mij„^([0,+oo[x[0,+cx)[). 

// both U and Z satisfy conditions (l)-(5) in the statement of Theorem \l.l\ then 
U{t,x) = Z{t,x) for a.e. {t,x) G [0, +oo[x [0, +cx)[. 

As a final remark, we point out tliat assuming tliat both the limits U and Z 
satisfy conditions (l)-(5) in the statement of Theorem 11.11 is reasonable. Indeed, 
concerning the classical viscous approximation (jl.lSp . these hypotheses are satisfied 
in the case when the global in time convergence has been proved, namely when 
B{U) = I, see Ancona and Bianchini [2] and also the analysis in Bianchini and 
Bressan [H Section 14]. As pointed out before, concerning the self-similar viscous 
approximation (|1.16p . convergence results have been established by Joseph and 
LeFloch [T71 [H] under suitable assumptions on the matrix B: their analysis shows 
that, under the same assumptions, condition (5b) holds. Conditions (l)-(4) in the 
statement of Theorem 11.11 are reasonable in view of the analysis in Tzavaras |31[ 
Sections 8-9] and Dafermos O Section 9.8]. 

The paper is organized as follows: the hypotheses exploited in the work are 
discussed in Section ILI) whereas in Section [2] we present some preliminary results 
that are used in the following section. The proof of Theorem 11.11 is established in 
Section |3] and Theorem 11.21 is an immediate consequence of Theorem 11.11 



1.1. Hypotheses. Here, we introduce the hypotheses of Theorem ll.il The first 
one is the standard assumption of strict hyperbolicity. 

Hypothesis 1. The flux F : M" — > M" is smooth and, for every U G W , the 
Jacobian matrix DF{U) has n real and distinct eigenvalues, 

\l{U)<\2{U)<---<\n{U). (1.17) 

Next, we also assume that the system of conservation laws is endowed with a 
strictly convex entropy. 

Hypothesis 2. There are functions ?/ : M" -^ M and q : W — > M such that 

Vj]{U) ■ DF{U) = Vq{U) (1.18) 

and D^r]{U) is strictly positive definite for every U G M". Here, " • " denotes the 
standard dot product. 

Last, we assume that the viscosity matrix B appearing in (|1.4p and (jl.Sp is 
dissipative. 

Hypothesis 3. For any given compact set H C K", there exists a constant an > 
such that, for every U G H , 

D^ri{U)v-B{U)v>aH\v\^ ViTgR". (1.19) 

Wc also quote Proposition 3 in Gisclon [11], since it is needed in Section |3l 

Lemma 1.3. // Hypotheses Q]^ hold, then, for every U G M", the number of 
eigenvalues of DF{U) with positive (resp. negative) real part is the same as the 
number of eigenvalues of B~^{U)DF{U) with positive (resp. negative) real part. 
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Now, we introduce the following standard terminology. The boundary in prob- 
lem (|1.1|) . ()1.7p is non characteristic if all the eigenvalues of the Jacobian matrix 
DF{U) are bounded away from zero, namely if there exists a constant c > such 
that 

Ai(C/) < • • ■ < A„_p(C/) < -c < < c < A„_j,+i((7) < • • • < \n{U) (1.20) 

for every U G R". It is actually sufficient to assume that (|1.20p holds for every 
U such that |[/ — J7o| < 3C(5, where C and 5 are the same constants as in the 
statement of Theorems 11.11 and 11.21 

Conversely, the boundary is characteristic if ()1.20|) is violated, namely if one 
eigenvalue of DF{U) can attain the value 0. If the boundary is characteristic, there 
arc positive constants A', c > satisfying 

Ai([/) < • • ■ < Xk-iiU) < -c < < c < Xk+iiU) <■■■< \n{U) (1.21) 

and 

|Afe(t/)|</« (1.22) 

for every value of U such that \U — Uq\ < 3C6, where C and 5 are the same constants 
as in statements of Theorems 11.11 and 11.21 

2. Preliminary results 

In this section, we collect some existing results that we exploit in the following 
section. In Subsection 12. II we give some results on invariant manifolds for ordinary 
differential equations and in Subsections 12.21 12.31 and 12.41 we consider self-similar 
solutions of systems of conservation laws. More precisely, in Subsection 12.21 we 
investigate the structure of a general self-similar, distributional solution having 
small total variation. In general, such a solution is not unique, so in Subsections l2.3l 
and 12.41 we describe how to construct a self-similar, distributional solution which 
has small total variation and which satisfies additional admissibility conditions. In 
particular, in Subsection 12. 31 we deal with the Riemann problem (|1.1|) . (|1.3[) . while 
in Subsection 12.41 we focus on the boundary Riemann problem (|l.ip . (|1.7p . 

2.1. Invariant manifolds for ordinary differential equations. In the follow- 
ing we rely on the notions of center, center-stable and uniformly stable manifold. 
For completeness, we present here the main properties that are exploited in the 
present paper. For an extended discussion, we refer the reader to the books by Ka- 
tok and Hasselblatt [^Hj and by Perko [17] and to the lecture notes by Bressan [7]. 
Let us consider the ordinary differential equation 

^ = G{V), (2.1) 

ay 

where V{y) G R*^ for every y and G : R"^ ^> R** is a smooth function which attains 
the zero value at some point. Without any loss of generality, we can assume that 
zero is an equilibrium point, namely G(0) = 0. 

Let us assume that at least one eigenvalue of the Jacobian matrix DG{0) has 
zero real part. Also, let us denote by V^ the subspace of R'^ spanned by vectors, 
the so-called generalized eigenvectors, associated with these eigenvalues and let us 
denote by Uc the dimension of V^. The Center Manifold Theorem ensures that there 
exists a so-called center manifold Ai''' C R**, which enjoys the following properties: 
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first, A4'^ is locally invariant for (|2.ip . namely, if Vb G M'*, then the solution of the 
Cauchy problem 

' dV/dy = G{V) 
1^(0) = Vo 

satisfies V{y) G Ai""' if l^l is small enough. Also, there exists a small enough ball, 
Br{Q), centered at 0, such that, if V{y) G Br{0) for every y, then V{y) G A^*^ for 
every y. Finally, M-'^ has dimension ric and it is tangent to V^ at 0. 

Let us denote by V the subspacc of W' spanned by the generalized eigenvectors 
of DG{0) associated with the eigenvalues with strictly negative real part. The 
Center-Stable Manifold Theorem (see e.g. [201 Chapter 6]) states that there exists 
a so-called center-stable manifold TM"^*, tangent to V (SV^ at the origin and locally 
invariant for (|2.ip . which satisfies the following properties: there exists a sufficiently 
small ball Br{0) such that, if V{y) is a solution of (|2.ip satisfying V{y) G 5^(0) for 
every y > 0, then V{y) G M'^^ for every y. 

Finally, let us consider the case that there exists a manifold £ which consists of 
equilibria of (|2.1I) . namely G{V) = for every V € £. Let E be the tangent space to 
£ at the origin. Then there exists a manifold A^^*, the so-called uniformly-stable 
manifold, which satisfies the following properties: first, A4|* is locally invariant 
for (|2.ip and it is tangent to V^ G) E at the origin. Also, there exists i?r(0), a small 
enough ball centered at 0, and a constant c > such that the following holds: if 
V{y) is a solution of (|2.ip satisfying 

lim \V{y) - Vle^y/^ = 

y— J--t-oo 

for some V G -8,(0) n £, then V{y) G AI^''' for every y. The uniformly-stable 
manifold can be viewed as a particular example of slaving manifold] see again |20j 
for a related discussion. 

2.2. Self-similar solutions of systems of conservation laws. Systems of con- 
servation laws are invariant under uniform stretching of both space and time vari- 
ables, and, thus, admit self-similar solutions in the form 

Uit,x) = Vi^) (2.2) 

where V{^) is a measurable function of one variable. The above form implies that 
self-similar solutions arc constant along straight-line rays emanating from the origin 
and, therefore, the study of such solutions is very important in the context of the 
Riemann problem, for which the initial data introduce a jump discontinuity at the 
origin. In this section, we consider self-similar solutions to Riemann problems and 
describe the structure of the solution. In particular, in this paper we consider initial- 
boundary value problems, so here we focus on the domain (t, x) g]0, +cxd[x]0, +oo[. 
Also, we assume that V has small total variation, i.e. 

TotVary < C6, (2.3) 

where < S << 1 and C > are as in the statement of Theorem 11.11 For a 
discussion on the main properties of functions in one space variable having bounded 
total variation, see Chapter 3.2 in the book by Ambrosio, Fusco and Pallara [1]. 

Having bounded total variation, V admits at most countably many discontinu- 
ities, say located at the points {jm}, m G N. The left and right limits of V, as ^ 
tends to a discontinuity point jm, exist and are finite and we adopt the standard 



UNIQUENESS CRITERION FOR THE BOUNDARY RIEMANN PROBLEM 9 

notation V^(jm) and V{j^). Also, the limits as ^ — >■ 0+ and £_ -^ +cx) exist and are 
finite and we denote them by 

lim V{C) = U and lim V{0 = Uq (2.4) 

respectively. Following Dafermos [U Section 9.1], we decompose ]0, +oo[ in disjoint 
sets as follows: 

]0,+00[=CU>VU U Jrn, (2.5) 

rnGN 

where C is the complement of the support of the measure dV/d^ and W is the 
(possibly empty) set of the points of continuity of V that are not in C. Since 
the measure dV/dS, vanishes on C, then the function V is constant on any interval 
included in C. From (|2.3p . we also have 

\ViO - Uo\ < C5 for every ^ e]0,+oo[. (2.6) 

In the remaining part of this section, we focus on the boundary characteristic 
case (|1.21[) - (|1.22[) . The non characteristic case (|1.20p can be treated in an analogous 
way. Let \i{V) be the i-th. eigenvalue of the Jacobian matrix DF{V) and by ri{V) 
a corresponding eigenvector with unit norm (the orientation is not important here) . 
Here, we introduce some notation that we need in the following. We define 

Aj+ := max XAV) + M5, A;- := min XAV)-M5 for i = fc + 1, . . . , n, 

\V-Uo\<CS \V-Uo\<C& 

Afc+ := max XdV) + M5. 

\V-Uq\<CS 

(2.7) 

Here, the constant M depends only on the flux F and the initial data Uq and its 
value is determined in the proof of Lemma 12.11 

Lemma 2.1. Let U{t,x) = V{x/t) be a self-similar solution of the system of con- 
servation laws 

dtU + d^F{U) =0, {t,x) g]0,+oo[x]0,+oo[ 

and assume that V satisfies p.3|) . // the system is strictly hyperbolic and the con- 
stant S in p.3|) is small enough, then the following hold. 

(1) There are constant states Ujl, . . . ,f/„_i G R", such that 

^^^>-[ Uo A„+<e<+oo ^'-^^ 

(2) At any point of discontinuity jm, ttt, e N, the Rankine-Hugoniot conditions 
hold, namely 

F{{V{j+)) - F{{j-,)) = j;„ [V{j+) - V{j-,)] for every m £ N. (2.9) 

(3) Let W be the same as in ()2.5|) . then, for every ^ G W, 

K{V{0)=^ (2.10) 

for some i = 1, . . . , n. IfVis also differentiate at ^, then dV/d^ is parallel 

tor,{ViO). 

(4) // ^ is in the (possibly empty) inner part of W and 

VA,(l-(O)-r,(F(O)^0, 
then V is differentiate at ^. 
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Proof. Properties (2), (3) and (4) follow straightforwardly from the analysis 
in [9l Section 9.1] and from Proposition 8.5 in Tzavaras [31]. Indeed, in the proof 
of that proposition one does not exploit that V is obtained as a limit of a viscous 
approximation. Regarding the structure (|2.8p . we assume that the constant S is so 
small that the intervals ]0, Xk + [, ]Afc+i — , Xk+i + [,•■• ]A„ — , A„ + [ are disjoint. One 
can first observe that from (|2.10p . it follows 

WC]0,Afe + [U]Afc+i-,Afc+i + [U ••■U]A„-,A„ + [. (2.11) 

Also, from the Rankinc-Hugoniot conditions (|2.9p . we get that 

' Df{9V{3 + ) + [1 - 9]vu-))de ■ [V{j + ) - V{j~)] ^ j„ [V{j+) - v{j-) 

and hence, jm is an eigenvalue of the matrix 

J DF{eVij + ) + [l-9]ViU)d9. 

This implies that the constant M in (|2.7p can be chosen in such a way that 
U Jm C]0,Afc + [U]Afc+i-,Afc+i + [U •■•U]A„-,A„ + [, 

mGN 

and, recalling (|2.1ip . Property (1) follows. ■ 

2.3. The solution of the Riemann problem. Let us consider the Riemann 
problem for a system of conservation laws, namely 

dtU + d^FiU)=0 

U+ x>0 (2.12) 



where {t,x) £ [0,+oo[xM and U{t,x) S M". Here, £/+ and U~ are two given 
constant states in R". In [22], Lax first introduced an admissibility condition on 
shocks, which is now named after him. Also, he constructed a distributional solution 
of (j2.12p satisfying this criterion. The analysis in [22] relied on the assumptions that 
|C/+ — U^ I is small and that each characteristic field is either genuinely nonlinear or 
linearly degenerate. The main point is the construction of the so-called i-wave fan 
curve through a given state C/"'", which in the following is denoted by Ti{si,U'^), 
Si being the variable parameterizing the curve, i = 1, . . . , n. The curve Ti{-, t/+) 
attains value in K" and satisfies the following property: given Si small enough, the 
Riemann problem 

' dtU + d^F{U)=0 

TT(n \-i U^ ^ > 

U{U,x)-<^ T,(s„[/+) a;<0 

admits a self-similar solution V{xlt) = U{x, t), which is either a rarefaction, a single 
contact discontinuity or a single shock satisfying the Lax admissibility condition. 
Also, let Xi{U) denote the i-th eigenvalue of DF{U), then there exists a small 
constant S > such that function V{£,) is identically equal to Ti{si,U'^) on the 
interval ^ g] — oo, Xi{U~^) — S] and is identically equal to U~^ on the interval ^ S 
[KiU+) + S,+oo[. 

In the general case, when U^ ^ Ti{si,U^), a solution of (|2.12p is obtained as 
follows: one defines the map 

X(5i, . . . , s„, U+) := ri(si, T2(52, Tal- • • , r„(s„, U+) ■ ■ ■ ))) (2.13) 
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and shows that it is locally invertible in a small enough neighborhood of (si, . . . , s„) = 
0. Hence, given U~ , by imposing U~ — x(si, . . . , s„, U~^), one uniquely determines 
the value of (si, . . . , s„). The corresponding solution of p.l2p is made by a finite 
number of rarefaction waves, contact discontinuities and shocks that satisfy the Lax 
admissibility condition. 

In [23l [24], Liu extended the analysis of Lax in [22] to very general systems 
and introduced the so-called Liu admissibility condition. Given C/"*" € M", the so- 
called i-Hugoniot locus is a curve Wi{s,U^), which is defined in a small enough 
neighborhood of U^ and comprises all the states such that the Rankine-Hugoniot 
conditions 

F{W,is, [/+)) - F([/+) = a,{s, U+) [W,is, U+) - U 

are satisfied for a certain real number ai{s, U'^) close to Xi{U~^). Given a value s, 
the shock between Wi{s, C/+) (on the left) and U~^ (on the right) is admissible in 
the sense of Liu if 

ai{s, [/"*") < (j{s, U~^) for every s between and s. (2-14) 

It should be mentioned that Liu admissibility condition can be viewed as an ex- 
tension of Lax admissibility condition. See Dafermos [S] Chapter 8.4] for a longer 
discussion. 

In U, Bianchini and Bressan constructed the solution of ()2.12p as the limit 
e — ^ 0"*" of the vanishing viscosity approximation, namely the limit of solutions W^ 
to the parabolic system 

" dtU' + d^F{U') ^ ed^^W 

U+ .T > (2.15) 



^ ' ' 1 [/ .T < 

Since in this work we exploit ideas and techniques from [l], we now go over the 
construction of the z-wave fan curve therein. If Si < 0, then one considers the 
following fixed point problem, defined on the interval [si , 0] : 

U^{t)^U++ I f,{U,iz),viz),a,iz))dz 

v{t) = fir) - conv/(T) (2.16) 

(7i{T) = — conv/(T). 

In the above problem, fi is the so-called generalized eigenvector defined in [U Sec- 
tion 4] by relying on the Center Manifold Theorem. Also, / is the generalized flux 
defined by 

/(t) = / Xi{Ui{z),v{z),a^{z))dz, 

where A^ is the generalized eigenvalue 

\i{Ui,v,ai) = DF{U)fi{Ui,v,ai) ■ fi{Ui,v,ai) . 

Finally, conv/ represents the convex envelope of /, given by 

conv/(r) = sup {(7(r) : g convex and g{z) < f{z) Vz G [si,0]}. (2-17) 

One can show that the fixed point problem (|2.16p admits a unique solution ([/^(t), v{t), (Ji (r)) 
and hence the i-th wave fan curve Ti{-,U^) is defined by setting Ti{si, t/+) := U{si). 
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If Si > 0, the construction is entirely similar, the only difference being that one has 
to take in (|2.16p the concave envelope of / instead of the convex one. 

The basic idea in the construction of the i-wave fan curve in [5] is the following: 
let us fix Si < and consider the solution of the fixed point problem (|2.16p . Assume 
that v{t) > on a given subintcrval ]a, /3[C [si, 0], and that v{a) — v{/3) ~ 0. Then, 
one can show that a solution of the Ricmann problem 

dtU + a,F(c/) - 

U{/3) x>0 




is a shock (or a contact discontinuity) which is also admissible in the sense of Liu. 
Conversely, assume that v{t) = on a given subintcrval ]a, b[C [sj, 0], then one can 
show that a solution of the Riemann problem 

d^F{U) = 

U{b) x>0 
U{a) X <0 

is cither a rarefaction wave or a contact discontinuity which is admissible in the 
sense of Liu. The solution is a contact discontinuity if VXi{U) ■ ri{U) = on ]a, h[. 
In general, the state U{si) is connected to C/+ by a sequence of rarefaction waves 
and of shocks (or contact discontinuities) satisfying the Liu admissibility condition. 
Once the wave fan curves Ti, . . . , T„ are defined, the argument in [4] works as in 
Lax [22]: one defines the function x as in p.l3p and shows that it is locally invert- 
ible, hence by setting C/~=x(si,...,Sn, U^) one uniquely determines (si, . . . , s„), 
provided |[/+ — U~\ is small enough. The limit e — !■ 0+ of (|2.15p is then given by 



U{t,x) 



In the previous expression, Ui and cti are obtained by considering the solution 
of the fixed point problem (|2.16p . It should be noted that in the case that each 
characteristic field is either genuine nonlinear or linearly degenerate, the wave fan 
curve constructed by Bianchini and Brcssan in [4] is the same as in Lax [22j . 

In [3], Bianchini extended the construction in 0] by considering more general 
approximations of the Riemann problem p.l2[) . In particular, he handled the 
viscous approximation 



u- 


X < cri(si) t 


U^{si) 


Cfi{0)t <X< (Ti+l(Sj+l)t 


Ur{r) 


X — cri(r)t i = 1, . . . , n 


U+ 


x>a„{0)t. 



dtU' + d^ 



F{U') ^ d^ B{W)d^W , (2.18) 



under some hypotheses on the matrix BiU). He also established a uniqueness 
result for (|2.12p . namely that there exists a unique solution U of ()2.12p which is 
obtained by patching together a sequence of rarefaction waves, shocks and contact 
discontinuities in such a way that the corresponding speed is increasing as one 
moves from the left to the right and in such a way that all shocks and contact 
discontinuities are admissible in the sense of Liu. In particular, this implies that the 
solution of ()2.12|) obtained by taking the limit of the viscous approximation (|2.18p 
does not change if one changes the viscosity matrix B{U). 
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Also, Dafermos [S] Chapter 9] provided a characterization of the hmit e ^' 0+ 
of the self-similar viscous approximation 



dtU' + d., F{U 



etdxxU^ 



^ y^^"^) \u- x<Q 

using a center manifold technique. In particular, he constructed the i-th wave fan 
curve for the self-similar approximation and showed that it is the same as the i-th 
wave fan curve Ti for the viscous approximation (|2.18p . See also the analysis in 
Tzavaras [30l [31] and Joseph and LeFloch [18] . 

2.4. The solution of the boundary Riemann problem. Consider the viscous 
approximation of a boundary Riemann problem 

dtU' + dx [F{U')\ = edx [B{W)dxU^] 

U^{0,x)^Ua, x>0 (2.19) 

U''lt,0)^UD, i>0, 

where Uq and Ud are two given constant states in R". For now, we focus on the 
case when the viscosity matrix B{U) is invertible for every U. 

The analysis of the boundary Riemann problem was initiated in Dubois and 
LeFloch [To]. One of the main difficulties in studying the limit of p.l9p is the 
presence of boundary layer phenomena. Convergence results have been proved in 
some special cases, see e.g. Gisclon [11] and Ancona and Bianchini [2]. 

The way boundary layers come into play is the following. Assume that C/^ con- 
verges as £ — > 0+ (in a suitable topology) to a limit function U having bounded total 
variation. One does not expect that U satisfies the boundary condition in (|2.19p . 
namely, in general, lim3,^o+ U{t,x) ^ Ud- Note that this limit is well defined since 
the solution U is assumed to be of bounded total variation with respect to the x 
variable. Let U denote the trace 

U= lim U(t,x), 

then one expects that the relation between U and Ud is the following: there exists 
a boundary layer W satisfying 

B{W)W' = F{W) - F{U) 

W{0) = Ud and Mmy^+^Wiy) ^ U. ^^■^''> 

Gisclon [11] and Gisclon and Serre [12] established, among other results, a detailed 
analysis of the boundary layers for a general viscous approximation when the vis- 
cosity matrix B is invertible. See also Joseph and LeFloch [TB] for the analysis of 
the boundary layers coming from other approximations. 

Concerning the analysis of the limit of (|2.19p . the case of identity matrix, B{U) = 
I, was considered by Ancona and Bianchini [5], who established convergence results 
and provided a characterization of the limit. The analysis of the limit was then 
extended in Bianchini and Spinolo [S] to the case of more general viscosity matrices, 
including the case that B{U) is singular, which is natural in view of most of the 
physical applications. For the convenience of the reader, we present the analysis 
in [5], since we use it in Section[3j In the case that the viscosity matrix B is singular, 
the initial-boundary value problem ()2.19|) may be ill-posed and hence one has to 
use a more complex formulation of the boundary condition. Here, for simplicity, we 
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discuss only the case of an invertible viscosity matrix. First, we focus on the case 
of a non-characteristic boundary, namely we assume that all the eigenvalues of the 
Jacobian matrix DF{U) arc bounded away from 0. The boundary characteristic 
case is considered afterwards. 

Let U denote as before the hyperbolic trace, then we expect that U and the 
Cauchy datum Uq are connected by a sequence of rarefaction waves, shocks, and 
contact discontinuities having positive speed and being admissible in the sense of 
Liu. Namely, 

U = Tn-p+l{Sn-p+l, T'n-p+2(Sn-p+2, Tn-p+ai' ' ' ' Tn{Sm Uq) ■■•))) (2-21) 

for some values s„_p_|_i, . . . , s„. In the previous expression, p denotes the number 
of positive eigenvalues of DF{U), while Ti{si, Uq) is the same i-wavc fan curve as in 
Bianchini [3] . We recall that the construction of Ti is presented in Section 12.31 To 
express the relation between U and Ud, we exploit (|2.20p and the Stable Manifold 
Theorem. This theorem ensures that there exists a manifold, the so-called stable 
manifold, having dimension equal to the number of eigenvalues of B~^(IJ)DF{U) 
with negative real part. Under some assumptions on B (see Lemma 11.31 in Sec- 
tion ll.ip . this number is exactly n—p. Moreover, all the functions W{y) satisfying 
B[W)W' = F{W) - F{U) and lim^^+oo W{y) = U lie on the stable manifold. 
Hence, if be denote by ip{si, . . . , s„_p, U) a function parameterizing the manifold, 
by combining the above argument with (|2.2ip . we get that 

Ud = ipisi, . . . ,Sn-p,Tn-p+i{- ■ ■ ,T„(s„,C/o) •••))))■ (2.22) 

Under some assumptions on the matrix B, the function 

(si, . . .,s„) 1-^ '/'(si, . ..,s„_p,r„_p+i(- ■ • ,T„(s„, C/o) ■ ■ ■)] 

is locally invertible and hence, by imposing p.22p . we determine uniquely the value 
of (si, . . . , s„). By arguing as in Section [2T3l the value U{t, x) of the limit of ()2.19p 
can be determined for almost every {t,x). 

Let us now consider the limit U of (|2.19p in the boundary characteristic case, 
namely we assume that the fc-th eigenvalue of the Jacobian matrix DFiU) can 
attain the value 0, i.e. Xk(U) ^ 0. The analysis in this case is much more delicate 
because there is not strict separation between the boundary layers on one side and 
the rarefaction waves, the shocks and the contact discontinuities with positive speed 
on the other side. 

For the boundary characteristic case, the construction in [5] works as follows. 
There exists a state, denoted by f/|, which is connected to Uq by rarefaction waves, 
shocks and contact discontinuities of the families k + I, . . . ,n: 

ul = Tk+i{sk+i,Tk+2{--- .Tn{sn,Uo) ■■■)), (2.23) 

for some Sk+i, . . . , s„. To complete the analysis, instead of the fc-wave fan curve 
Tfc, one employs the characteristic wave fan curve 1k{sk, Ul), which is constructed 
as follows. If Sk < 0, we consider the fixed point problem 

t/fc(r)=C/«+ / h{Ukiz),Vkiz),a{z))dz 







VkiT) ~ /(t) — monconv/(r) (2.24) 

air) ~ — -— monconv/(r), 
ddr 



UNIQUENESS CRITERION FOR THE BOUNDARY RIEMANN PROBLEM 15 

where t e [sfc,0]. Here, rk is the same as in p.l6p and / denotes again a "gen- 
eralized flux" . Finally, monconv/ represents the monotone convex envelope of /, 
namely 

monconv/(T) = sup{(7(t) : g convex, non decreasing, g{z) < f{z) Vz £ [sfc,0]}. 

(2.25) 
By relying on the Contraction Map Theorem, one can show that there exists a 
miique continuous solution {Ukir), ^^(t), (t(t)) of p.24|) which is confined in a small 
enough neighborhood of (f/^, 0, Afe(t/|)). One eventually defines the characteristic 
wave fan curve by setting Tfe(sfe,/7|) :— U{sk)- If Sfe > 0, the construction is 
entirely similar, the only difference being that in (|2.24p one has to take the monotone 
concave envelope of / instead of the monotone convex one. 

The idea behind the construction sketched above is the following. Let us fix 
Sfc < and consider the solution of (|2.24[) . defined on the interval [sk, 0] and set 

s :— maxJT G [sfc, 0] : (t{t) = 0} 

and 

s := supjr e [sfc, 0] s.t. (T(r) = and tife(r) > 0}. 
Clearly, s < s. Let us focus on the most interesting case by assuming Sk < s_ < 
s < 0. Then U{s) (on the left) is connected to U{0) = Ul (on the right) by a 
sequence of rarefaction waves and shocks (or contact discontinuities) satisfying the 
Liu admissibility condition and having strictly positive speed. Also, U{s) (on the 
left) is connected to U{s) (on the right) by a sequence of contact discontinuities 
that are admissible in the sense of Liu and have speed 0. Moreover, by construction, 
Vk{T) > and ct{t) = for t £ [sfc,s[. Hence, the Cauchy problem 

dy ^ 1 

dT vkir) (2.26) 

y{s) = 

defines a change of variables: t i-^ y, [s/t, s[— )• [0, -|-oo[. One can show that Wk{y) '■= 
Uk{T{y)) is a solution of the equation B{W)W' = F{W) - F{Uk{s)) satisfying 
Wfc(O) = %k{sk, ul) and Wuvy^+o. Wk{y) = Uk{s). 

To complete the construction, one defines a manifold, transversal to the vector 
1k{sk, ul) — ul, having dimension k — 1 and satisfying the following property: if 
VFs(O) lies on this manifold, namely Ws(0) = i^(si, . . . , Sk-i) for some parameteri- 
zation Lo of the manifold, then the solution Ws of the equation 

B{W. + ul)[Ws + Ulr]' = [F{W.. + ul)]' (2.27) 

decays exponentially fast to the equilibrium point as y — > -|-cxd. Moreover, for 
fixed (si, . . . , Sk-i) and Sk, one singles out a function Wp, such that by setting 

Wiy) = Wsiy) + Wkiy) + Wj,iy), (2.28) 

one obtains a solution W of system [_B(W^)W^']' — [F{W)]' satisfying limy_j.+oo W{y) = 
Uk{s), where Uk is the solution of (|2.24p . In (|2.28|) . Wk{y) is the above-mentioned 
function obtained by considering the change of variables (|2.26[) and Wg (y) is a so- 
lution of (|2.27|) satisfying W^s(O) = uj{si, . . . , Sk-i). One can determine Wp in such 
a way that Wp is small (in the sense specified in [5]) with respect to Ws and Wk, 
namely it can be regarded as a perturbation term. By construction, W(jj), and 
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hence the value M^(0), depends on the vector (si, . . . , Sk-i), on the number Sk and 
on the state C/|. One sets 

M^(0):=^(si,...Sfc,C/f) (2.29) 

and then, by recalhng (|2.23p . one defines the map 

(si,...,s„,[/o) >^ tp(si,...,Sk-i,Sk,Tk-i{--- ,T„(s„, t/o))j, 

which, under suitable assumptions on the matrix B, is locally invertiblc. Hence, by 
imposing 

Ud = 1p{si,...Sn,Uo), 

one uniquely determines (si, . . . , s„) and hence the limit of (j2.19p . This completes 
the construction in [S] for the boundary characteristic case. 

Finally, we underline that in [ITl [18] , Joseph and LcFloch had previously con- 
structed a solution of the boundary Riemann problem (|l.ip . (|1.7p . by taking the 
limit e — > 0"'^ of the self-similar viscous approximation 

dtU' + d^ [F{U')\ = etd^ [B{U')d^U'] 
C/(0, x) = Uo and U{t, 0)^Ud, 

and they established the convergence under suitable hypotheses on B. Also, they 
provided a precise description of the limit, in both the non-characteristic and the 
characteristic case. In the boundary characteristic case, most of the analysis relied 
on the assumption that the ^-characteristic field is genuinely nonlinear. 

3. Proof of Theorem 11.11 

The existence part of the theorem is proved in [5j Section 3] . More precisely, by 
applying Theorems 3.1 and 3.2 in |5], we get the following. There exist constants 
C and 5i , Si sufficiently small, such that if 

\Uo-Ud\<Si, (3.1) 

then there exist a function V satisfying conditions (l)-(5) in the statement of 
Theorem 1 1 . 1 1 with S = 6i. Here, we focus on the issue of uniqueness. We distinguish 
between two cases: if the boundary is characteristic, namely if one eigenvalue of 
the Jacobian matrix DF{U) can attain the value 0, then the proof relies on two 
preliminary results. Lemmas 13.11 and 13.21 which are introduced in Sections 13.11 
and 13 . 21 respcctivelv. In the case when the boundary is non characteristic, the proof 
of Theorem 11.11 is actually much easier and it follows from an argument in Joseph 
and LcFloch [17]. Both cases are discussed in Section [373] 



3.1. A decomposition result. In this subsection, the principal result is a decom- 
position of the boundary layers and its construction is given in Lemma 13.11 

Throughout this subsection and the following one, we focus on the boundary 
characteristic case, i.e. we assume that 

Ai(C/) < • • • < Afc_i(t/) < -c < < c < \k+i{U) <■■■< Xn{U) (3.2) 

and 

\Xk{U)\<KSi, (3.3) 

for some constants K > and the same constant di as in p. II) . 
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To begin with, we introduce some notations: we set 

A{U, U') = DF{U) - \b{U)\ ' 

and consider boundary layers, namely functions W satisfying 

W" = B-'^{W)A(W, W')W' , ^. 

limj,_)._i_oo W{y) exists and it is finite. ^ ' ' 

The above ordinary differential equation can be rewritten as a first order system 

W =p 

p' = B-\W)[A{W,p)p-ap] (3.5) 

a' = 0. 

The reasons why we introduce the "fake unknown" a would become clear in Sec- 
tions 13.21 and 13.31 

For any given state C/" G M", we can linearize p.Sp about the equilibrium point 
(U^,0,Xk{U^)). By applying the Center Manifold Theorem, we infer that there 
exists a center manifold A4''' having dimension n + 2. For some more details of the 
construction, we refer the reader to Step 1 in the proof of the lemma below. Note 
that the center manifold is not unique, but we can arbitrarily fix one. 

Lemma 3.1. Let conditions (|1.17p - p.l9|) in Section \l.l\ be satisfied and let C/" he 
a state in M" satisfying jC/" — C/qJ < C5i. Then there exist constants Ci, C2 and 5, 
(5 < <5i, such that the following holds. Assume that (|3.2p and p.3p hold, and that 
W is a solution of (j3.4p satisfying 

\W{y) -U^< CS, for every y e [0, +cx)[ 

limy_,+oo W'{y) = (3.6) 

\W'{y)\<C6 Vye[0,+oo[. 

Then W can be decomposed as 

W{y) = Uk{y) + U,{y) + Up{y) for every y € [0, +«)[, (3.7) 

where 

(1) Uk{y) lies on the center manifold Ai'^ for every y £ [0, +oo[ and satisfies 
Vimy^+ooUl{y) =0. 

(2) Us satisfies 

\Us{y)\ < CiSe-'-y^^ for every y £ [0, +cx)[ (3.8) 

and B{Us + U»)U'J{y) - A(C/, + [/«, U'^)U'^. 

(3) the perturbation term Up is small with respect to the previous ones, more 
precisely 

\Upiy)\ < 025^6-"^/^ for every y e [0, +oo[. (3.9) 

The constant c in p.8p and p.9p is the same as in 



Proof. We exploit the construction in [2] and [5] and we rely on the notions of 
center, center-stable and uniformly stable manifold. For an overview, see Section r2.1l 
and the references therein. The proof is established in several steps. 
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Step 1 We linearize system (|3.5p about the equilibrium point ([/", 0, Afc ([/")) and 
we obtain the (2n + 1) x (2?! + 1) matrix 

/„ 

i3-i(;7«)[A([/«,0)-Afc([/«)/„] 

0* d' 

where 0* is the zero row vector in R" and and /„ denote respectively the null and 
the identity matrices with dimension n x n. By relying on Lemma ll.3[ we get that 
the matrix 

i3-i(;7«)[A([/«,0)-Afc(t/«)/„ 

has the same number of eigenvalues with non-positive real part as the matrix 
[A{U^,0) — \k{U^)In)- Consider the center-stable space, that is, the subspace 
of K^"+^ generated by the generalized eigenvectors corresponding to eigenvalues 
with non positive real part. Because of p.2p and p.3p . the center-stable space 
has dimension n + 1 + k and therefore, by the Center-Stable Theorem, there ex- 
ists a {n + 1 + fc)-dimcnsional center-stable manifold Ai'^'^ containing any solu- 
tion of (|3.5p that, for y G [0, -l-oo[, is confined in a small enough neighborhood of 
([/",0,Afc([/*')). Let us fix a function W satisfying ^^^ and ([321): if S is small 
enough, then (W{y),p{y) = W'{y),0) lies on M'^^. 

Step 2 In the following, we use the basis composed by the eigenvectors n,. . . ,rk of 
A{U^ ,Q) in M". It should be noted that these eigenvectors are linearly independent 
because, by assumption, A has n real and distinct eigenvalues. By arguing as in [5j 
pages 45-47], we get that W'{y) admits the representation 

W'{y) = R,s[w{y),V,s{y),())v,s{y), (3.10) 

where Vcs{y) is a function taking values in M.^ and -Res is a function taking values 
in the space of n x fc matrices and satisfying the following: Rcs{U\Q^\k{U'^)) is 
the matrix whose columns are ri • • • r^ . 

Step 3 We now employ the same "diagonalization procedure" as in [5J pages 47-51] 
and wc write 

VUy) = Rs{w{y),VM,Qys{y) +fk{w{y),Vk{y).Q)vk{y), (3.11) 

where Vs and Vk take values in M.^~^ and R, respectively . The functions f^ and Rg 
take values in M'^ and in the space of fc x (fc — 1) matrices respectively. Also, from 
the construction in [5] it follows that, for every W and Vk, the point 



(VF, Rcs{W, 0, Vk, Q)fk{W, Vk,Q)vk, 



lies on the center manifold M'^ defined just before the statement of Lemma [3?T 
Moreover, for every W and Vg, the point 



{W, Rcs{W, F„ 0, {))Rs{W, Vs,0)Vs,Q 



belongs to the uniformly stable manifold Al|'* of system p.5|) with respect to the 
manifold of equilibria 



f = /(W^,0,ct), W ^W\ae 



We refer the reader to Section 12.11 for the notion of uniformly stable manifold. Fi- 
nally, from the construction in [5] follows that fk (/7*, 0, Afe([/'')) ^ ru and Rs (t/», 0, \k{U^)) 
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is the matrix whose columns are ri ■ ■ ■ Vk-i- Hence, if W{y) is fixed, the map 
{Vs, Vk) t-5- Vcs is locaUy invertible. In view of the above analysis, we consider W{y) 
to be the solution of p.4p and we first obtain the function Vcs (y) from (|3.10p and 
then, using the local invertibility of the map, we get the functions Vs{y) and Vk{y) 
satisfying (|3.1ip . 

Step 4. By relying again on the "diagonalization procedure" in jSJ pages 47-51], 
we get that {W{y),Vk{y),Vs{y)) satisfy 

W = RUW,Vkfk + RsVs,0)fkVk + Rcs{W,Vkfk + RsVs,0)RsVs 
v'^=MW,Vk,Vs,0)vk (3.12) 

V^ = AiW,Vk,Vs,0)Vs, 

where (f>k is a real- valued function satisfying <f)k{U^,0,Q,Xk{U^)) = 0, while A attains 
values in the space of (fc — 1) x (fc — 1) matrices and A{U^, 0, 0, Xk{U^j) is a diagonal 
matrix having entries Ai(C/'') . . . Afe_i(C/''). 

Step 5. We now illustrate our strategy to complete the proof of the lemma, the 
details are provided in the following step. 
We write 

Vs{y) = Vs{y) + Vp{y) and Vk{y) = Mv) + Mv)^ (3-13) 

where (T4,Wfe) arc the "principal components" and {Vp,Vp) the "perturbations". 
We impose that the components Us and Uk in (|3.7p satisfy 

Uk = Rcs{Uk,Vkfk,0)fk{Uk,Vk,0)vk , . 

i'^^MUk,ikJ,0)vk ^ ■ ^ 



and 

U'^ = RcsiUs + UKRsVs,0)RsiUs + C/«, Vs,0)Vs 
V^ = A{Us + 0,O,Vs,O)Vs 



(3.15) 



respectively. 

Hence, by exploiting dSj]) and ([3?T3| and by plugging ([3?T4| and (|3J5l) in ([3?T2l) 
we obtain the equation for {Up,Vp,Vp). Then, by imposing that both Up and Us 
converge to as y — s- +oo, we get a fixed point problem for (Us, Up,Vp, Vp), which 
by applying the Contraction Map Theorem is shown to admit a unique solution (in 
a suitable metric space). 

We then get Vs = Vs ~Vp, Vk ~ Vk — Vp and Uk — W — Up — Us- by construction, 
Uk and Us satisfy the properties described in the statement of the lemma. 
Step 6 We now provide the details of the argument sketched in Step 5. One has 
to keep in mind that now W{y), Vs{y) and Vk{y) are given functions, and that 
from (j3.13p one gets the identities Vs = Vs — Vp and Vk = Vk — Vp- Hence, by 
integrating the first line in p.l5p one gets that the equation for Us{y) 

Us{y)^ I Rcs{Us + uKRs{Vs-Vp){z),0)Rs{Us + uKVs-Vp,Q){z)\Vs{z)+ 

— Vp{z) dz. 

(3.16) 
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Combining equations p.l2|) - p.l5|) . we arrive at the equation for Up{y) 



Upiy) 



-(-CX3 



R^cs{W{z),Vkfk{z) + R,V,{z),0)fk{z)+ 

- Rcs{W{z) - Us{z) - Up[z),Vkfkk[z),Q)fkk{z) Vk{z)dz+ 

y 

Rcs{W{z) - Us{z) - Up{z),Vkfkkiz),0)fkkVp{z)dz~ 

Rcs{W{z),Vkfk{z) + RsVs{z),0}Rs{z)+ 



+ 00 

y 



+00 



(3.17) 



Rcs{Us{z) + UKVsRss,0)Rss{z) V,{z)dz + 



+ 00 



Rcs{Usiz) + uKVsRs..,0}R,siz)Vpiz)dz. 



To simplify the exposition, we empfoy the notations fk ~ rk{W, Vk,0), fkk = fkiW— 
Us - Up, ik, 0), Rs = Rs{W, F„ 0) and i?,, = i?,([/, + 0, 14, 0). Similarly, we get 
the equation for Vp 



y 

+00 



iJpiy) ^ / 4>k{W - Us ~ Up,Vk ~ Vp,0,0){z)vp{z)dz + 

J +00 

[(f>k{W, Vk, Vs,0) ~(t)k{W~Us- Up, Vk ~ Vp, 0, 0)(z)] Vk{z)dz , 

(3.18) 
and the equation for Vp 

AiW,Vk, Vs,0) - A{Us + t/«,0, Vs - Vp,0)\ Vsiz)dz + 



Vpiy) = I e 

'0 



y _ 



My- 



Here, we have set 



AiW,Vk,Vs,0)-A Vp{z)dz. 



A = A(c/«,0,0,Afc([/«)), 



(3.19) 



(3.20) 



namely A is a diagonal matrix whose entries \i{U^), . . . , \k-i{U^) have all strictly 
negative real part. 

We now introduce the sets 

Xs := {Us e C"([0,+oo[;]R") : \Us{y)\ < C^Se-'^l^ for every y} 
X^l := {Up e C°([0, +^[; K") : \Up{v)\ < C26^e'^y/^ for every y} 
Xl :- {Vp e C°([0,+oo[;R^-i) : \Vs{v)\ < C^b^e-'^'^ for every y} 
XI -^ {vp £ C°([0, +oo[) : \vpiy)\ < C^S^er^y/^ for every y}, 

where c is the same as in (|1.2ip and Ci, C2, C3 and C4 are suitable constants to 
be determined. We equip each of the previous spaces with the norm 

||F|U:= sup {\F{y)\e^y/^}, 
ae[o,+oo[ 
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SO that Xs, Xp , Xp and X^ are all complete metric spaces. We also equip the 
product space Xs x X^ x X^ x X^ with the norm 

\\iUs,Up,Vp,Vp)\\ := \\U,\U + ||C/p|U + \\Vp\U + \\vp\U. 

On XsXX^xX^x X^ we define the map T = (Ti, T2, Tg, r4) by setting Ti, Ta, T3 
and T4 equal to the right hand sides of p.l6p . p.l7p . p.l9p and p.lSp respectively. 
We also recall that from p.lOp we deduce that 1 14,5(2/)! ^ ^(1)<5 f^^ every y. 
Hence, by applying the Local Invertibility Theorem to the map (T4,i'fc) i-> Vcs 
defined as in p. lip , we obtain the functions Vg and v^ satisfying 

\Vs{y)l \vkiy)\ < C5S , Vj/e[0,+oo[, (3.21) 

for a suitable constant C5 > 0. 

We now introduce the following estimate, which is proved in STEP 7: there 
exists a constant Cg > such that 

1^5(^)1 < Cee-^y^^S , for every y > 0. (3.22) 

By relying on p.22p and p.2ip . one can then show that, if S is small enough, 
then one can choose the constants Ci , C2 , C3 and C4 in such a way that T attains 
values in Xg x X^ x Xp x X^ and is a strict contraction. The computations are 
quite standard and are similar to (but easier than) those performed to prove Lemma 
3.9 in [5], so we skip them. 

Finally, by applying the Contraction Map Theorem we infer that the systems 
obtained by combining equations (|3.16p . p.l7p . p.l9p and p. 181) admit a unique 
solution belonging to X^ x Xp x Xp x X'^. 

STEP 7. We conclude by proving p.22p . Let us consider the set 

X:={yeC°([0,+oo[,R'=-i) : \\V\\s<Ced}, 
where the value of the constant Cg will be determined in the following and 

\\V\\s= sup {e^y^'\Viy)\}. 

ye[0, + oo[ 

It turns out that X, equipped with the distance induced by the norm || ■ ||s, is a 
closed metric space. Also, if the constant Cg is big enough, then the map defined 

by 

G{V){y) = e^«K(0) + T i^^^-'^ \A{W{z),Vk{z), V{z),0) - a] V{z)dz 
Jo L J 

is a strict contraction from X to itself. To prove that G is a contraction, one 
exploits bounds (|3.6p and (|3.21l) . The fixed point of G is the solution of the ODE 
V' = A{W, Vk, V, Q)V which attains the value K(0) at ?; = 0, hence from ((3T2)) we 
have V{y) = Vs{y) for every y. From the definition of X it follows that Vg satisfies 
the bound (g^. ■ 

Remark 3.1. As the above proof shows, the heuristic idea underlying Lemma IXTj 
is that we can "invert" the construction in [2] (or in [5l pages 45-68]). Indeed, 
in [2] one first finds (Uk,Vk) and (Us,Vs) by solving (|3.14p and p.lSp respectively 
and then obtains (C/, Ufe,!^) by constructing the perturbation term {Up,Vp,Vp). 
Conversely, in Lemma \3A] we are given {U,Vk,Vs) and we want to obtain {Uk,Vk), 
{Us,Vs) and {Up,Vp,Vp). 
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3.2. Uniqueness of the characteristic wave fan curve. In this subsection, 
we establish the uniqueness of the characteristic wave fan curve for the boundary 
characteristic system. It should be noted that the construction of the characteristic 
wave fan curve is presented in Section 12.31 

Let U{t, x) ~ V{x/t) be a self-similar solution of the conservation law (|1.12[) hav- 
ing small enough total variation. Then V has the structure described in Section [221 
and in particular there exist constant states Ul, . . . ,t/„_i £ M", such that 

^^(t\^fuf Xi+ <£, < Xi+i-, i^k,...,n-l . ^ 

^^^^^l uo A„+<e<+oo, ^■^■^■^> 

where A^-l- and A^— are defined in (j2.7p . Also, in the following we use the same nota- 
tion as in (|2.4p . we denote by U the right limit of F at 0, namely U ~ limj^o+ ^(Oj 
and we recall the constants Si and C in the beginning of Section [3] 

Lemma 3.2. Let the functions F and B satisfy conditions (J1.17p - (|1.19[) . Then 
there exists a sufficiently small constant 6, < S < 6i, such that, if (|3.2p and (|3.3p 
hold and ifU{t,x) — V{x/t) is a self-similar solution of the conservation law (|1.12p 
such that 

(1) TotVarV<CS, 

(2) all the shocks and the contact discontinuities of U are admissible in the 
sense of Liu, 

then the following holds. Let Ul and U he as in (|3.23p and (|2.4p respectively and 
assume that Ub is a state satisfying 

(3) there exists a value U_ such that the following two properties hold: 

(a) IJ_ (on the left) and U (on the right) are connected by a shock (or a con- 
tact discontinuity) having speed and satisfying the Liu admissibility 
condition; 

(b) there exists a solution of (|3.14p . entirely contained in a ball centered 
at {Ul,0) and having radius CS. Such a solution satisfies Uk{0) = Ub, 
limj^^+oo Uk{y) = U and linxy^+oo ^/(.(y) = 0. 

Then Ub ~ 1k{s,Ul) for some s, where Tfc is the characteristic wave fan curve 
defined in ()2.24p . 

Proof. The proof exploits many ideas and techniques of Bianchini from [3J 
Theorem 3.2]. An overview of his analysis can be found in Section [^^ The main 
difference is that we also need to treat the presence of the boundary layers. To deal 
with this, we employ the monotone convex envelopes of the generalized flux. 

The proof is divided into several steps: 
Step 1 Fix the orientation of the eigenvector rk{U) in such a way that 

[U^Ul]-rk{Ul)<Q, 

where as before the symbol "•" denotes the standard dot product in M". By relying 
on the analysis in [3] , we infer that there exists a value s < such that the solution 
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of the fixed point problem 

U{t) = Ul+ f fk{Uiz),viz),aiz))dz 

JQ 

«(r)=/(T)-conv[,,o]/(T) (3.24) 

1 dconv[g,o]/(T) 

"^") = d Tr 

satisfies U{s) — U. Here, convjjo] denotes the convex envelope of the function / 
computed on the interval [s, 0] and the function / is defined by 



fir) 



^k{U{z),v{z),0,a{z))+da{z] 



dz + A, 



(3.25) 



where A is an arbitrary constant. Indeed, the solution of (|3.24p does not depend 
on the choice of A. In p.25p . (j)^ is the same function as in p.l2[) and the constant 
d is given by 



d:=- 



bu{U,vu,Vs,a) 



da 



U=Ul,Vk=0,Vs=0,cr=\k{Ul) 



(3.26) 



and it is strictly positive by construction (see [51 formulas (3.32)-(3.33)]). 
Step 2 By combining the assumption (3a) in the statement of the lemma with the 
analysis in [3l Section 3], we deduce that there are two cases, the first one being 
that there exists s<s such that the fixed point problem 



U{t) = U+ h{U{z),v{z),a{z))dz 



«('^) = fiT) - conv[^,s]/(r) 
1 dconv[^_5]/(T) 



(3.27) 



a{r) 



dT 



defined on the interval [s, s] with 

fir)= f k(C/(z),t^(z),0,a(z))+dCT(z) 



dz 



(3.28) 



has a solution {U, v, a) that satisfies U{s) = U_. 

The second case is that there exists a s > s such that the solution of the fixed 
point problem 



U{t)^U+ Pk{U{z),v{z),a{z))dz 



v{t) = 1{t) - conC[5^s]£(T) 

f dconcr,5,5i^(r) 
cr(r) - 



(3.29) 



d dr 

satisfies U{s) ~ U_. Note that the fixed point problem is now defined on the interval 
[s, s] . In the previous expression, the function £ is defined as in (j3.28p , namely 



lir) 



(j)k{U{z),v(z),Q,(j{z)) + d(j{z) 



dz. 



(3.30) 
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Figure 1. The case we rule out in Step 3 




Also, eonC[s sj£ denotes the concave envelope of i on [s, s], namely conC[s^sj^ := 
— conv[j s](— ^). In Step 3, this second possibility is ruled out. 
Step 3 The goal in this step is to show that the second case presented above cannot 
occur. To prove this, we proceed by contradiction and exploit Lemma 3.1 in |2]. 

Assume that p.29p defines a zero speed shock (or contact discontinuity) which 
is admissible in the sense of Liu and which connects IJ on the right to IJ_ on the 
left. In other words, assume that the solution (t/, w,(t) of (j3.29|) satisfies lJ{s) = C/, 
conC[s^sj^ = and 

1{t) < for every r e [s,s]. (3.31) 

Indeed, the fact that concjj s]£(r) = is equivalent to the admissibility of the shock. 
Let us now fix ft, < nnn{|s|, \s — s\} (see Figure [1]) and consider the map 



(3.32) 



T,[U,v,a]iT)^U+ fk{U{z),viz),aiz))dz 
T2[U,v,a]{T) = £(r) - conC[g,g+,,]£(T) 

rr^^rr .f X 1 dCOnC[g^g+fe] £(t) 

Ti[U,v,a]{T) = - — , 

defined on the interval [s, s + h]. Also, let us define the norm 

\\{U,v,a)\\h ■- C,\\U\\c»([s.s+h]) +C,\\v\\co([s^s+h]) + h\\L^([s,s+h]), (3.33) 

where C is a small enough positive constant to be determined. Then, we get that, if 
C is small enough, the map T = (Ti, T2, T3) defined in p.32[) is a strict contraction 
from a suitable metric space into itself with respect to the norm || • ||/i. In the 
following, we denote by 7'' = ([/'', w'', cr'') the unique fixed point of T and by i^ 
the quantity 

^''(T)=£([/^^;'^a'')(T) . 
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By relying on the same argument as in the proof of Lemma 3.1 in [3], we get that 

a^{T) < for every t e[s,s + h]. (3.34) 

Since p.30p imphcs that l{s) — 0, then conC[g^s+h]^'^{s) = 0. Hence inequahty p.34p 
imphes that 

e'^ir) < conC[5 g+,,]^''(T) < for every T€[s,s + h]. (3.35) 

Let us now return to p.24p : since (j{t) > for t g]s, 0], then by arguing as before 
we get 

< conv[5^o]/(^) < fir) 
for every r e] s, 0] . Hence p.35p imphes 

i'^ir) < fir) for every r G ]s, s + h]. (3.36) 

Here and in the fohowing, we denote by / the restriction to the interval [s, s + h] 
of the same function / as in p.24[) . Next, we estimate ||/ — ^^\\c°{[s.s+h])- First, 
we point out that, for every r € [s,s + h], we have 

f{T)-&{T)= f \mU,v,<j) - MU\v\<j'')]{T)dT + f d{a-a'')dT 



< 0{l)h\\U - W^WcOilsrs+h]) + Oil)h\\v ~ v''\\cO[s.s+h] + Clk - <T''\\LHls,s+h]y 

(3.37) 

Here and in the sequel, we denote by 0(1) any constant depending only on the 
flux F and on the state C/^.. Also, we exploit the estimate \d4>k/da + d\ < 0{1)S 
and, by taking S sufficiently small, we can assume \d(j>k/d(J + d\ < C^- By choosing 
/i < (^^ in (|3.37p . relying on the definition (|3.33p and using the fact that the map 
T = (Ti,T2,T3) defined in p.32p is Lipschitz continuous with constant 1/2 and 
that 7'' is the fixed point of T, we eventually get 

fir) - i\r) < 0(1)C^||7 - I'h < Y^h T(7)lk, (3-38) 

where 7 = (C/, v, a) is the restriction of the solution of (|3.24p to the interval [s, s + h]. 
To estimate ||/— /''||c''([s.s+W)i it remains to provide a bound on II7— T(7)||/i. Byre- 
calling that the solution of p.24p satisfies C/(s) = U and taking into account p.32p . 
we deduce 

[/ - Ti (7) = on [s, s + h] (3.39) 

and 

|w(t) - 72(7)(t)| = |conC[j,s+ft]/(r) - conv^j^oj/MI- 
By relying on the chain of inequalities 

max / > conC[s.s+h]fir) > fir) > conv[j.o]/(r) > > ^''(r) Vr e]s, s + h], 

(3.40) 
we obtain 

\v{t) - T2(7)(t)| < conC[s^j+,j]/(r) - conv[g^o]/(T) < max / 

Tels,s+h] /g ^^. 

< Il/-^''l!c0([s,s+/i])- 



Finally, let us consider the quantity 



h-T3{-f)\\Lmg.S+h]) 



s+h 



conv[5^o]/ conC[g_5+,i]/ 



dr dr 



dr . (3.42) 
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Since conv^gQ^ f /dr — coiiC[g^s+h]f /dr is non-decreasing, the set of zeroes is an in- 
terval, which can be empty or reduced to a point. Let us fix s* such that 



conv[s^o]f/dT - conC[s^s+h]f/dT < in [s, 



and 



conv[5_o] f/d-T - conC[g^5+^] //dr > in [s*,s + h]. 
Using p.42p . we estimate 



lk-T3(7)|[ii([j_5+,,]) 



conC[s ^s+h]f/dT - conv[j,o]//(^T 

s+h ^ 



dT~ 



= 2 



< 2 



conv[s^o]f/dT ~ conC[s^g+h]f/dT 
conc[,5,5+h]/-conv[5.o]/ (s*)+ 



dr 



(3.43) 



^5,0]/ ~ conC[g^g+h]f 



is + h)< 



conc[,5_5+,i]/-conv[5^o]/ («*) 



< 



< 2||/ - £'\\c°{ls,s+h])- 

To get the last two inequalities, we exploit (|3.40l) and p.4fp respectively. 
Combining (|339)) . ([QT]) and ((03)) with (|338)) . we arrive at 

ll/(r) -^''(r)||co([,.,+;,]) < 0(f)CII/(r) -^"(r)||co([,,,+;,]). 

By taking C sufficiently small, the above inequality implies ||/(t)— ^''(T)||po(rj j^/ji^ = 
0. Since this contradicts (|3.36p . we conclude that (|3.29p cannot hold. 
Step 4 From here and on, we investigate the first case stated in Step 2, namely the 
fixed point problem (|3?27| - (|3?28| . In this step, we "glue together" (|3^ and ([S^Tf]) 
by finding a fixed point problem defined on the interval [s, s] . 

First, we point out that the function /(r) defined by (|3.25p for r G [s, 0] and 
by p.28p for r G [s.,s[ is continuous on [s, 0], because in p.25p the constant A is 
chosen in such a way that f{s) = 0. We define the function g on [s, 0] by setting 

{conv[g.o]/(r) if r e [s, 0] 
(3.44) 
conv[3^j]/(r) ifTe[s,s[. 

Since we assume that U_ and U are connected by a shock (or a contact discontinuity) 
with speed zero, then dconv^g^f/dr is actually on [s,s]. Also, by assumption 
dconv[5 0]//'^''' is strictly positive on [s, 0]. It then turns out that the function g is 
convex. In addition, by construction ^(t) < /(r) in [s, 0], hence 

5(r) < conv[s o]/(''') for every r 6 [s, 0]. (3.45) 

Conversely, the restriction of convj^ o]/(''') to the interval [s, 0] is a convex function 
less than or equal to /, hence 

conv[s_o]/(T) < conv[g^o]/('r) for every t e [s,0]. 

By applying the same argument to the interval [s, s] , we deduce that 

g{T) > conv[^_ol/('^) for every r G [s,0] (3.46) 
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and therefore, by combining inequality p.45p with p.46p . we arrive at the identity 

g(r) == convi^,o]fiT) for every r e [s, 0]. 

Wc finally deduce that the solution of the fixed point problem (defined on the 
interval [s, 0]) 

' U{t) = uI + j h{U{z),v{z),a{z))dz 

v{T)=f{T)-conv^,o]f{T) (3.47) 

l<iconv[^.o]/(T) 

-"(") = d Tr 

satisfies U{s) = U and U{s) ~ U_. 

Step 5 We now show that, under the hypotheses of the lemma, one can use the 
monotone convex envelope in (j3.47p instead of the convex one. The monotone con- 
vex envelope is defined in (|2.25p . Let C£ ([s, 0]) denote the space of the functions 
/ defined on [s, 0] such that / is continuously difFerentiable functions and /' is 
Lipschitz continuous with Lipschitz constant bounded by L. If / G C^' ([s, 0]), 
then conv[s_o]/ G C!£ ([s, 0]). This follows from more general results discussed by 
Griewank and Rabier in [Ij. Also, if / S C^' ([s, 0]), then 

r/ m\ f COnv/(fs,0l)(T) if T > To fn An\ 

monconv/(s,0 ) = <^ r}\ Ji) \ -c J (3.48) 

•'^-' ^' \ conv/([s,0])(To) ifr<To, ^ ' 

where 

To :— min s t £ [s, 0] : -^conv/([s, 0]) > 
[_ dr 

If the derivative of the function conv/([s, 0]) is strictly negative in [s, 0], then one 
sets To = s and then, relation p.48|) is still valid. The proof of p.48|) can be found, 

for example, in [51 Page 361. Since in p.47p a — conv/([s, 01) > 0. then by 

a dr ' 

exploiting (|3.48p . we infer that instead of the fixed point problem (|3.47p . we can 

consider the following one: 

U{t) = uI+ f rk{U{z),viz),aiz))dz 





vir) = /(t) - monconv[s^o]/('r) (3-49) 

_ 1 (imonconv[s,o]/(T) 
d dr 

which is still defined on the interval [s, 0] . 

Step 6 We now exploit assumption (3b) in the statement of the lemma. Namely, 
there exists a solution ([/, Vk) of the system 

U' = Rcs{U,Vkfk,0)fk{U,Vk,0,0)vk ,„,„^ 

which satisfies [/(O) ~ Ub, limy_i.+oo U{y) = U_ and limy^.+oo U'{y) = 0. RccaU 
that Res and (j)k are constructed in the proof of Lemma [33] 
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li U_ — Ub, then the conclusion of the lemma follows from the previous steps. 
So in the following we assume U_ ^ Ub- Since any point of the form (L'^, 0) is an 
equilibrium for system (j3.50p . then Vk{y) 7^ for every y. 

Let us first consider the case when Vk{y) > for every y. The change of variables 



dr 
dy 



Vk{y) 



maps the interval [0,+oo[ into some interval [s^s[ and {U,Vk) satisfies 

U{t) = U+ I rk{U,Vk,6,0)dz 



In the previous expression, we employ the equalities ffc = Rcs^k and Imiy^^oo Vkiy) = 0. 
The last condition follows from the fact that U'{y) converges to and \fk \ is bounded 
away from 0. By setting 

/(t)= / MUiz),Vk{z)AO)dz, 

J s_ 

we have Vfc(r) ~ /(r) > for every r S [s,s[. Since /(O) = 0, then = 
nionconv[s,s[/ and hence we can rewrite the previous system as 



U{t) = U+ rk{U{z),v{z),a{z))dz 



«(t) = /(t) - monconv[s^s[/(T) 



(3.51) 



a{r) 



1 rfmonconvfj, ^[/(t) 
d dr 







By arguing as in Step 4, we can then "glue together" (|3.49p and p.5ip and obtain 
that the solution of the fixed point problem 



U{t)^uI+ h{U{z),v{z),u{z))dz 



w(t) = /(t) - monconv[s_o]/(T) 

_ 1 rfmonconv[g^o]/(T) 
d dr 

for r G [s, 0] satisfies U{s) = Ub- Following a similar argument as in Step 3, we 
can rule out the possibility that Vk{y) < in p.50p . The proof of the lemma is 
complete. ■ 



3.3. Conclusion of the proof. Here, we complete the proof of the theorem. In 
Subsection l3.3.1l the boundary characteristic case is treated and in Subsection l3.3.21 
the non-characteristic. 
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3.3.1. Boundary characteristic case. By relying on Lemma 12. ![ we infer that the 
self-similar function V satisfies 

V( f:\_iul Aj+<C<A,+i-, i = k,...,n-l , . 

^^^)-{ ul = Uo A„+<^<+oo. ^^-^^^ 

By relying on assumptions (3) and (4) in the statement of Theorem 11.11 we have 
that, for any i = k,...,n — l,Uf and U^^i are connected by a sequence of rarefac- 
tions and shocks (or contact discontinuities) that are admissible in the sense of Liu. 
Hence, by exploiting Theorem 3.1 in j3|, we get that, for every i = k, . . . ,n — 1, 

uI^TMM+i) 

for some Si small. In the previous expression, Ti denotes the same i-wave fan curve 
as in formula p.lSp . Hence, 

Ul = Tfc+i (.Sfc+i, rfc+2(. . . , T„(s„, Uo)... )) (3.53) 

for a suitable vector (s^+i, . . . s„). 

We now exploit assumption (5) in the statement of Theorem 11.11 and we apply 
Lemma [XT] to the function W . Let Uk{y) be the same function as in the statement 
of Lemma 13.11 then by setting Uk{0) = Ub and applying Lemma 13.21 we obtain 
Ub = Tk{Ul, Sk) for some small s^. By comparing the construction in [5l Section 
3.2.3] with the proof of Lemma l3.ll one infers that Ud ~ i^isi, ■ ■ ■ ,Sk,Ul), for 
a suitable vector (si,...,Sfc), where ip is the same function as in ()2.29p . Hence, 
relation p.53p leads to 

Ud = iIj(si, . . . ,Sk,Tk+i{sk+i,Tk+i{. . . ,Tn{sn,Uo) . ..))) = 0(si, . . . , s„, C/q)- 

The analysis in J5l page 67-68] ensures that the map (j) is locally invertible, hence 
given Ud and C/q such that \Ud — Uo\ is small enough, the value of (si, . . . , s„) is 
uniquely determined by the above relation. 

In conclusion, if [/ is a solution satisfying all the assumptions of Theorem 11.11 
then the value U{t, x) is uniquely determined for a.e. {t, x) G [0, -|-oo[x [0, -|-oo[ and 
it can be obtained as follows. Set as in p.52|) C/| ~ Uq and for any i ~ k, . . .n — 1 
define inductively Uf as Ti^i{si+i,Uf^^). Also, let {Ui,Vi,(Ji) be the fixed point 
of (I2.16|) (if Si < 0, otherwise one takes the concave envelope instead of the convex 
one). Also, assume for now that Sk < and denote by {Uk,Vk,o'k) the fixed point 
of (I2.24P and by s the value 

s = maxJT G [sfe, 0] : o-/c(t) = 0}. 
Then, 

if x/t = ak{T) tG[s,0] 
ri(f'r)=^ ^i a <7i{0) < x/t < at+i{si+i) foTi = k,...,n-l 

^' ' ^ '' '-^ if.T/t = aj(r) iovj=k + l,...,n 
\i x/t > cr„(0). 

If Sfe > 0, one needs to take in (|2.24p . the monotone concave envelope and define 
the value s as 

s = min{r G [0, Sk] : (Jkir) = 0} . 
The proof for the boundary characteristic case is thus complete. 
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3.3.2. Non characteristic boundary case. We now focus on the case when the bound- 
ary is non characteristic, namely we assume that (|1.20p holds. 

By arguing as in Scction r3.3.1l and relying on assumptions (l)-(4) in the statement 
of the theorem, wc infer that [7, the trace of U on the t axis, satisfies 

U = T„_p+ifs„_p+i,T„_p+2(. . . ,T„(s„,C/o) • ■ •)] (3.54) 

for a suitable vector {s^-p+i, ■ ■ ■ Sn). 

Let us now exploit the assumption (5) in the statement of the theorem: since the 
boundary is non characteristic, then / is locally invertiblc and from the relations 
fit)) = /([/) and \U ~U\ < CS, wc obtain U ^ U_. Hence, system (|1.14p becomes 

B{W)W' = F{W) - F{U) 

W{0) = Ud limj,^+oo Wiy) = U 

and by exploiting again the fact that DF{U) is invertible, we get that Ud lies on 
the stable manifold of the above system about the equilibrium point [7, namely 

[/£) = V3(si,...,s„_p,C7j (3.55) 

for a suitable vector (si, . . . , Sn-p) and for the same function tp as in (J2.22I) . 
By combining p.54p with (|3.55p . we arrive at 

Ud = tpl Sl, ■ • ■ , S„_p,T„_p+i( S„_p+i,T'r,_p+i(. . . ,T„{Sn, Uq) ■ ■ ■)) I . 

Once Ud and Uq are given, this relation uniquely determines (si, . . . , s„) and, thus, 
the function U satisfying assumptions (l)-(5) in the statement of the theorem. D 
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